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SUMMARY 

This  memorandum  is  concerned  with  the  motion  of  a  long,  flexible  wire 
in  orbit.  The  wire  is  approximated  by  n  point  masses  joined  by  massless 
rigid  rods. 

Under  the  assumption  of  an  inverse  square  gravitational  field,  the 
general  equations  of  motion  of  the  n  masses  are  developed.  Then,  with  all 
external  forces  except  gravity  neglected,  and  with  the  motion  of  the  wire 
restricted  to  the  orbital  plane,  two  specific  cases  are  examined.  In  the 
first  case  the  wire  spins  about  its  center  of  mass,  after  having  been 
initially  extended  in  some  fashion.  The  second  case  examines  the  motion 
of  an  extended  wire  which,  initially,  is  approximately  aligned  with  the 
vertical. 

The  analysis  indicates  that,  in  the  absence  of  internal  dissipative 
forces,  the  straight-line  configuration  of  the  spinning  wire  is  not  stable, 
since  a  parametric  resonance  condition  exists.  In  the  second  case  the 
wire  will  oscillate  about  the  vertical  if  the  initial  disturbances  are 


small. 
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vector  position  of  the  center  of  mass  with  respect  to  an  earth- 
centered.  inertial  reference  frame 

vector  position  of  the  i-th  particle  with  respect  to  the  system 
center  of  mass 

external  gravitational  force  acting  on  i-th  particle 
internal  system  forces  acting  on  the  i-th  particle 
mass  of  the  i-th  particle 
total  mass  of  system  of  particles 

coordinates  located  at  center  of  mass,  rotating  with  respect  to  the 
local  vertical 

angle  between  local  vertical  and  the  r)  axis  (see  Fig.  2) 
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I.  INTRODUCTION 

The  use  of  satellites  for  various  communication  purposes  has  been  under 
study  at  HAM)  for  some  time.  One  proposed  type  of  satellite  consists  of  a 
long  wire  in  orbit.  The  wire,  which  serves  as  an  antenna,  would  be  kept  in 
tension  by  spinning  it  about  its  center  of  mass.  This  device  would  provide 
intermittent  communication  between  a  sending  and  a  receiving  station. 

A  wire  is  attractive  for  this  application  since,  being  flexible,  a  long 
antenna  may  be  contained  in  a  relatively  contact  package  for  the  ascent 
phase. 

The  relationship  between  the  orbital  altitude,  the  spin  rate,  the 
frequency  of  the  transmitted  radiation  and  the  characteristics  of  the  antenna 
have  been  examined  in  Ref.  1,  while  in  Ref.  2  the  mechanical  design 
requirements  of  a  rotating  antenna  have  been  considered. 

System  studies  indicate  that  an  interesting  design  is  a  wire  approxi¬ 
mately  2000  ft  long,  orbiting  at  an  altitude  of  from  1500  to  2000  miles. 

For  these  conditions,  the  wire  should  not  deviate  from  a  straight  line  by 
more  than  x/8,  where  1  is  the  wavelength  of  the  incident  radiation/1^  For 
typical  values  of  X,  the  permissible  transverse  deflection  of  the  cable  is 
of  the  order  of  a  few  centimeters!^  Various  design  considerations  lead 

to  a  spin  rate,  with  respect  to  the  local  vertical,  of  approximately  one 

(l  2) 

radian  per  second.  9  This  research  memorandum  is  concerned  with  the 
stability  of  a  spinning  wire  when  subjected  to  alternating  tension  and  com¬ 
pression  forces  due  to  the  gravitational  gradient. 

After  the  packaged  wire  has  been  placed  into  orbit,  the  wire  must  be 
unwound  or  extended  in  some  fashion  and  then  accelerated  to  the  desired  spin 
rate.  Once  the  wire  is  spinning  in  the  orbital  plane,  it  may  oscillate  due 
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to  the  influence  of  the  gravitational  gradient  and  other  small  perturbing 
forces.  Thus  the  question  arises,  -whether  the  straight-line  configuration 
of  the  antenna  is  stable. 

In  this  research  memorandum,  the  spin-up  phase  is  not  considered.  It 
is  assumed  that  initially  the  wire  is  spinning  in  the  orbital  plane,  and 
that  the  transverse  deflections  are  very  small  compared  to  the  length  of  the 
antenna.  Since  a  wire  is  a  continuous,  flexible  body,  its  transverse  motion 
can  be  expressed  in  terms  of  partial  differential  equations.  However,  the 
gross  characteristics  of  a  flexible  wire  may  be  approximated  by  a  number  of 
point  masses  joined  by  rigid  rods  without  mass.  It  is  hoped  that  by  examin¬ 
ing  a  simple  model  of  this  sort,  any  obvious  stability  problems  will  be 
revealed.  Furthermore,  insight  should  be  gained  to  the  mathematical  and 
physical  approximations  that  can  be  made  in  an  analysis  of  a  continuous  wire. 
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II.  DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

In  Fig.  1,  the  variables  used  to  define  the  motion  of  the  center  of 
mass  are  indicated,  while  Fig.  2  shows  the  coordinates  used  to  define  the 
positions  of  the  n  point  masses.  Only  the  planar  case  will  be  considered. 


a 


Fig.  I  —  Center  of  mass  coordinates 
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Fig. 2 — Particle  coordinates 


The  acceleration  of  the  i-th  mass  with  respect  to  an  inertial  frame  is 


El. 

1 


+ 


F  .  +  F  . 
gi  si 


(1) 


where  F  .  is  the  external  gravitational  force  acting  on  the  i-th  particle 
and  F  .  represents  all  of  the  forces  internal  to  the  system  which  act  on 

S  X 

the  i-th  particle.  From  the  definition  of  the  center  of  mass 

n 

X  mi  pi  ■  0  (2) 

1 


Noting  that  ^  Fg^  must  be  zero,  Eq.  (l)  may  be  written  as  follows: 


M(r  )T  =  >  F  . 

o'l  gi 


(3a) 


n 

-  -r  L 


F  .  +  F  .  +  F  . 
gi  gi  si 


(3b) 
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Equation  (5a)  defines  the  motion  of  the  center  of  mass  of  the  wire,  while 
Eq.  (3b)  defines  the  motion  of  the  i-th  particle  with  respect  to  the  center 
of  mass  as  observed  from  an  inertial  reference  frame.  When  the  motion  of 
the  i-th  particle  is  seen  from  the  rotating  t]-£  coordinate  system,  Eq.  (jb) 
has  the  following  form: 

m. 
l 

=  -  ~  Vf.+f.+f. 

M  gi  gi  si 

1 


(4) 


+  2cd 


ms 


x  (p.  )  + 

'ms 


(co  )  x  p.  + 
ms  ms  l 


CD  x 
ms 


(CD  X 

ms 


pi} 


where  the  subscript,  ms,  indicates  that  vector  differentiation  is  with 
respect  to  the  rotating  system,  and  to  is  the  angular  velocity  of  the  tj-£ 
system  with  respect  to  the  inertial  frame.  Writing  Eq.  (4)  in  terms  of  its 
scalar  components  yields 


m. 

i 


i±  -  2ri.  (a  +  G)  -  T)i(a  +  e)  -  (a  +  e)2 


+  F 


+  F 

gC1  s^i 


(5a) 


m. 

i 


n 

.  ....  .  .  p  "I  V”' 

Tii  +  2^  (a  +  ©)  +  £.(a  +  e)-n.  (a+er  =  •  iT 


gill 


(5b) 


+  F  .  +  F  . 
grii  sr|i 


Thus  for  n  particles  there  are  2n  equations  identical  in  form  to  Eq.  (5). 
However,  all  of  the  2n  equations  are  not  independent,  because  of  center  of 
mass  considerations,  and  because  the  total  length,  l,  of  the  ’wire’  constrains 
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where  i  runs  from  1  to  n-1. 

Thus  the  internal  force  components  acting  on  the  i-th  particle  may  be 
written  as  follows: 


F  . 

ST|1 

-  FTi,i-l  C0S 

Pi-l,i 

-  FTi,i+l  COS 

(8a) 

F  y. 

=  FTi,i-l  SlD 

^i-1,1 

FTi,i+l  Sln  Pi,i+1 

(8b) 

where  F^  i+1  is  a  tension  force,  and  =  FT,i+l,i’  Tiie  tension 

forces  are  unknown  forces  of  constraint  which  arise  due  to  the  accelerations 
of  the  connected  particles. 

Introducing  Eqs.  (7)  and  (8)  into  Eqs.  (5a)  and  (5b)  yields: 
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•  P  "I  m  •  v — ' 

!i  li  -  (“  +  «)  -  \  (a  +  e)  -  5l  (a  +  e)  ‘  -  ^  L  V  +  Fgti 


+  F  f  '  h  )  .  F  f  _klii±i 


n 

i±  \  +  2^  (a  +  9)  +  q  (a  +  9)  -  ri,  (a  +  9)2  =  -  ~  ^ 

- 


F  .  +  F 
gr(i  gT]X 


+  F  (  ’ll~1  '  ^ 

Ti>1-1  V  W 


11,1+1  \  h.m 


Finally,  the  external  gravitational  force  vhich  acts  on  the  i-th  particle 


may  he  determined  from  consideration  of  Fig.  o. 


Fig.  3 — Resolution  of  gravitational  force  into 
the  ?7-£  coordinate  system 
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F 

gCi 


n.  g  r„  r 
x  °o  E  o 

3 

r . 
i 


sin  a  + 


(10a) 


F  . 
gT)i 


m.  g  r  r 
i  o  E  o 


cos  a  + 


J 


(10b) 


The  distance  from  the  center  of  the  earth  to  the  i-th  particle  may  be  found 
from  the  law  of  cosines : 


2  2  2 

r.  =  r  +  p.  +  2r 
X  o  'l  o 


T) .  cos  a  +  t.  sin  a 
x 


(11) 


Thus 


-3  ~  -o 
r .  =  r 

1  o 


3  ,  1i  h 

1  -  ^  (2  —  cos  a  +  2  —  sxn  a. 
2  r  r 

o  o 


(12) 


In  the  binomial  expansion  only  first-order  terms  have  been  retained. 
Introducing  Eqs.  (10a)  and  (10b)  and  (12)  into  Eqs.  (9a)  and  (9b)  yields: 


;n2 


2  n.  (a  +  ©)  -  rii  (a  +  9)  -  i±  (a  +  0) 


(l3a ) 


So  rE 


2 

(  -  3r)i  sin  a  cos  a  -  sin  a  + 


+  FTi,i-l  f  ^i-1  _  FTi,i+l  /  ^i  ~  ^i+1  \ 

mi  V  ^i-l,i  /  mi  V  i,i+l  / 


\  +  2^  (a  +  ©)  +  q  (a  +  6)  -  r)i  (a  +  0) 


go  rE 


(13b) 


- —  (  -  3^  cos  a  -  3^  sin  oc  cos  a  +  tj^  ) 


+  f  r‘i-l,i  ^i\  FTi L i+1  j  ^i  ^i+l 

mi  \  \-l,i  /  mi  \  ^iji+l 


The  influence  of  higher  order  terms  is  considered  in  Ref.  3. 
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Thus  for  n  particles  interconnected  by  (n-l)  massless,  rigid  rods 
there  are  2n  equations  of  the  form  of  Eqs.  (13a)  and  (13b)  and  (3n-l) 
unknowns.  The  remaining  (n-l)  equations  are  provided  by  Eqs.  (7a)  and  (7b). 

1  i,i+l  =  ^i  "  ^i+l^  +  ^i  "  ^i+1^  (14) 

i  =  1, . . . ,  n-l 

The  manner  in  which  the  t)  -  £  axes  rotate  with  respect  to  the  local 
vertical  has  not  yet  been  considered.  One  reasonable  approach  is  to  select 
a  as  a  function  of  time  based  on  the  motion  of  the  equivalent  rigid  body. 
Rather  than  prescribing  the  rotational  motion  of  the  *1  -  £  axes  it  is  also 
possible  to  constrain  one  of  the  tj  '  s  or  f^'s  by  making  a  an  unknown.  The 
latter  procedure  has  been  followed  in  the  subsequent  development  of  the 
equations  of  motion. 

It  is  assumed  that  the  number  of  mass  particles,  n,  is  odd.  If,  in 
the  undeflected  state,  the  'cable'  exhibits  mirror  symmetry  with  respect 
to  a  plane  perpendicular  to  its  center,  then  there  will  be  a  particle  at 
that  center.  By  the  use  of  Eq.  (14)  all  of  the  Vs,  with  the  exception  of 
one,  may  be  expressed  in  terms  of  the  inter-particle  lengths  and  the  Vs* 

The  remaining  r)  will  be  taken  to  be  that  of  the  center  particle  and  its 
value  will  be  constrained  to  be  zero.  Thus  for  n  particles  there  are  2n 
equations,  n  unknown  £'s>  (n-l)  unknown  tension  forces,  and  a.  However, 
the  constraint  imposed  upon  Tln+1  does  not  result  in  a  unique  set  of  T)  -  £ 
axes  (see  Fig.  4).  For  the  symmetric  distribution  of  particles  illustrated 
in  Fig.  (4),  either  the  tj  -  £  or  the  tj'  -  £ '  axes  would  satisfy  the  imposed 
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Fig. 4 — Possible  77  —  £  axes 


constraints.  For  the  undeflected  cable  to  always  coincide  with  the  t)  axis, 
the  following  condition  must  be  satisfied: 


n-1  n 


Equation  (15)  insures  that  the  tj-£  axes  are  such  that  center  of  mass  of 
one  half  of  the  cable  has  a  mirror  image. 

Finally,  in  view  of  the  magnitude  of  the  permissible  transverse 
deflection  as  discussed  in  Section  I,  only  terms  containing  the  first  power 
of  ^  will  be  retained.  As  a  consequence, 

i,i+l  ^i  ^i+l  (16) 

Equation  (16)  together  with  Eq.  (6a)  indicates  that  all  of  the  derivatives 
of  the  q's  must  be  zero. 

Taking  into  consideration  Eqs.  (6a),  (6b),  (15)  and  (16)  it  is  now 
possible  to  obtain  from  Eq.  (13a)  an  expression  for  the  angular  motion  of 
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the  q  -  £  axes.  This  is  done  by  summing  first  over  the  top  half  of  the 
cable  and  then  subtracting  from  that  result  the  summation  over  the  bottom 
half. 


where  the  index  j  is  identical  with  — —  .  If  the  "cable"  were  rigid,  all 
of  the  £'s  would  be  identically  zero,  and  Eq.  (17)  then  reduces  to  the 
form  that  is  usually  considered  when  studying  gravitational  gradient 
effects. 

Following  the  same  slimming  procedure  with  respect  to  Eq.  (13b)  yields 


If  we  only  consider  the  motion  of  the  center  particle,  Eq.  (13b)  becomes 
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2  £  (a  +  e)  +  (a  +  e) 


5m,  g„  rT 


j  &o  E  £ .  sin  a  cos  a 

5  J 

r 

o 

+  "  FTj,j+l 


(19) 


However,  the  components  of  the  tension  force  acting  on  the  j-th  particle 
must  be  equal  in  magnitude  to  the  summations  of  the  tension  components 
acting  on  the  remaining  particles. 


J-l 

I 


1  L 


F  -  F 

rTi, i-1  Ti,i+1 


=  -  F, 


(20a) 


I 

J+l 


FTi,i-l  ‘  FTi,i+l 


'TJ,j+l 


(20b) 


J-ir 


I 


FTi,i-l  ^  ^i-1  "  ^i ^  "  FTi,i+l  ^i  "  ^i+l) 


'  FTj,j-l  (^j-l  '  ^J) 


I 

J+l1- 


FTi,i-l  ^i-1  ’  ^  "  FTi,i+l  ^i  “  ^i+1^ 


FTJ,j+l  ^0  "  ^j+1* 


(20c) 


(20d) 


Thus  the  tension  forces  which  act  on  the  j-th  particle  may  be  determined 
in  terms  of  the  variables  a  and  by  consideration  of  Eqs.  (l8),  (19)> 


(20a)  and  (20b). 
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FTj,J-l  “j 


(a  +  0)  +  i  (a  +  6) 


/  ^‘1 

^  I  ”i  ”i )  (a  +  e)‘ 


(21a) 


5g^  r 


g„  rx 


°  E  m.,  sin  a  cos  a  -  ■  °  — E  (1-5  cos2  a)  f  X  mi 


2r 


5  2  SJ 


c 


TJ 


,J+i  = (a  +  6)2  ( I  “i  ’’i) 


3g  r| 

+ - —  m.  t  sin  a  cos  a 

i  -*  -i 


2r 


5  2 


(21b) 


-  m 


Cj  (a  +  e)  +  i  5j  (a  +  0) 


go  rE 


J  r: 


2  /  j-1 

(1-5  cos2  a)  (  Y,  mi 


Utilizing  Eqs.  (20c),  (20d),  (21a),  and  (21b),  Eq.  (17)  may  be  written  as 
follows : 


5g, 


a  +  0  + 


sin  a  cos  a 


-  —L_  (t 
21  ^j-l 
P 


-  W 


(a  +  0)2  - 


gn  P 

-2-—  (l-  5  cos2  a) 


(22) 


where  1  is  the  inter -particle  distance,  assumed  to  be  a  constant  for  all  of 
the  particles. 

Introducing  Eqs.  (20a)  and  (20b)  and  (21a)  and  (21b)  into  Eq.  (15s) 
yields  the  following  transverse  equation  of  motion  for  the  center  particle. 


(a  +  0)2  -  -2-^  (1 


-  5  cos^  a) 


-  (a  +  0)‘ 
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go  rE  M  ,  .  2  , 

+  ■  -----  (1-3  sin  a) 


j-1 

Y  m.  rj. 

r  i  i 


1  m . 
P  J 


(a  +  8)‘ 


(23) 


°  E  (1-5  cos2  a) 


+  Sj+P 


We  are  now  able  to  write  the  remaining  transverse  equations  of  motion  in 
terms  of  functions  of  a  and  the  C’s.  The  tension  force  acting  on  the  j-1 
particle  due  to  the  j-2  particle  is,  from  Eq.  (13b) 


FTj-l,j-2  PTj-l,J  +  mj-l 


2£  (a  +  e)  +  (a  +  e) 


(24) 


•  o  ^  r 

A\2  .  O  E 


-  T]^  (a  +  6)  +  (rij.i  -  3^^  sin  a  cos  a  -  3^_1  cos 


a) 


However,  Fm,  .  is  identical  with  F_,  .  ,  of  Eq.  (21a).  Substituting 
Eqs.  (21a),  (22)  and  (24)  into  Eq.  (13a)  yields  the  following  expression 
for  the  transverse  motion  of  the  j-1  particle: 


^.1-1  +  ^  j-1 


(a  +  6)  - 


2 

-  _E  (1-3  cos2  a) 


J-1 

X  m. 

_ * 

m .  ..  1 

J-1  P 


% 


3  T) 


j-1 


21 


P  -I 


-  (a  +  e)2 


+ 


2 

r 
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(1-3  sin2 


(a  +  e)2 


(25) 
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2 

&  r  P 

-2 — -  (1-5  cos2  a) 


J-l 

Z 

1 


mi  \ 


Vi  S 


^j-2  +  ^ 


Vl 

21 


(Cj+l  +  2^J- 


The  procedure  that  has  been  followed  in  obtaining  Eq.  (25)  may  be 
applied  in  turn  to  the  j-2  particle,  etc.  Because  of  the  symmetry  that 
exists,  the  equation  of  motion  of  any  one  of  the  particles  comprising  the 
lower  half  of  the  'wire*  may  be  determined  from  the  corresponding  upper-half 
equation  by  substitution  of  indices,  and  noting  that 


nj-l  “  ^j+l'  6tC* 


(26) 


From  center  of  mass  considerations 


m  £  = 

n  ^n 


1  > 
’  5  mJ 


-L  u 

-  2  mj 


m£  +  ...  +  ?j.j. 


m ...  t  m  ,  £  , 

j+l  bj+l  n-1  bn-l 


(27a) 

(27b) 


Thus  after  the  elimination  of  the  tension  force  expressions,  there  are  (n+l) 
equations  to  be  solved.  These  are  the  (n-2)£  equations,  the  a  equation  and 
the  two  equations  of  the  center  of  mass  motion.  Of  these  equations,  (n-l) 


are  coupled. 
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III.  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 


From  Fig.  1  and  Eqs.  (5a)  and  (12),  the  orbital  equations  of  motion 


are : 


i2 

r  -  r  6  = 

o  o 


go  rE 


(28a) 


d  (r2  e)  =  0 


dt  o 


(28b) 


Thus  the  rotational  and  vibrational  degrees  of  freedom  do  not  influence 

* 

the  center-of-mass  motion.  It  will  be  assumed  that  the  orbit  is  circular. 


and  that  6  equals  cn  ,  a  constant. 

The  equations  to  be  solved  are: 


••  2 

a  +  >c  sin  a  cos  a  =  gj —  -  CJ+1) 


(29a) 


+  ^j*k 


‘  5a  +  6i  <■ 


2Q  +  m,  .  t) 


j-k  'j-kv 

B.  1  ' 


j“k  P 


=  S. 


Q  +  m.  ,  tj 


-  r  w  /  j-k  j-k  \  f 
21p  UJ-1  W  1  lp  '  5j*(k-l) 


+  iiiio^i) 


m .  .  1 

j-k  p 


(29b) 


*  _  5 

This  is  due  to  the  fact  that  in  the  binomial  expansion  of  r.  ,  only 

first-order  terms  have  been  retained.  (See  Ref.  5.) 
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where 


k  —  0)  I  •  •  • ;  (j-2) 


■  2  2  2 
s  =  (a  +  co  )  -  a>  (1-3  cos  a) 


s2  =  (a  +  CDC)2  -  co2  (1-3  sin2  a) 


Q  = 


j-l 

I 


j-k 


mi^i  “ 


j-l 


m.  Tj , 

1  1 


j-k 

Z  m. n.  =  0  if  j-k  >  j-l 

.  ,  i  1 

J-l 


In  Eq.  (29b),  negative  values  of  the  index  k  correspond  to  particles 
located  in  the  plus  tj  half  of  the  tj  -  £  plane,  and  the  converse  is  true 
for  positive  values  of  the  index. 

An  examination  of  Eqs.  (29&)  and  (29b)  reveals  a  set  of  (n-l)  coupled, 
non-linear  homogeneous  equations.  It  is  immediately  obvious  that,  for 
values  of  a  on  the  order  of  cdc,  the  stability  of  the  system  of  particles 
is  questionable.  For  this  condition  the  alternating  tension  and  compression 
forces  due  to  the  gravitational  gradient  are  predominant. 

Two  limiting  cases  may  be  analyzed  by  approximate  methods.  The  first 
case  is  that  in  which  the  spin  rate  of  the  wire  is  very  large  compared  to 
o>c.  In  the  second  case  the  wire  is  almost  aligned  with  the  direction  of 
the  local  vertical,  and  the  angular  rate  with  respect  to  the  vertical  is 
very  small. 

Case  I 


We  will  assume  that  a  is  almost  a  constant. 
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a  =  Q  +  A  d 


(30) 


where 


n  =  d(o)  »  A  a 

The  parameters  s^  and  s^  may  be  written  as  follows: 

p  n  p 

s.  ~  n  +  2  a  AQ  +  2  n  CO  +  2co  Aa  +  ~  of  (1  +  COS  2a)  (31a) 

1  c  c  2  c 


s0  =  fl2+2QAa+2fio  +  2co  Aa  +  \  cu2  (l  -  cos  2a)  (311) 

2  cede 


For  orbit  altitudes  ranging  between  300  and  2000  miles,  cu  is  of  the 

order  of  10~^  radians  per  second.  With  2  approximately  equal  to  one 

2 

radian  per  second,  s^  and  s^  are  almost  equal  to  0  .  Thus  it  might  be 

expected  that  the  form  of  the  solution  of  Eq.  (29b)  is  primarily  oscilla¬ 
tory  in  nature.  Since,  in  the  derivation  of  the  equations  of  motion, 
the  £'s  have  been  assumed  to  be  small,  products  such  as  £  Aa  will  be 
neglected  as  higher  order  terms.  With  this  assumption  Eqs.  (29a)  and 
(29b)  may  be  approximated  as  follows: 


*  ^2 

Aa  +  ^  cd  sin  2a  = 
2  c 


^j-i  "  jj+i; 
21 


o  2 

+  22  cdc  +  a)c  (l  +  cos  2a) 


(32a) 
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^j+k  + 


ft2  +  2  ft  CD  +  ^CD^(l  -  COS  2a) 
c  2  c 


n£  +  2  (1  0)  +  |  CD2  (l  +  cos  2a)  1  ^  - ^~k  ^"k 

c  2  c 


’]( 


m .  .  1 

J-k  P 


} 


ft  +  2  ft  CD  +■ 
c 


■|  cd2  (l  +  cos  2a)J 


iifis  (t  -  t  ) 

21  Vbj-1  bj+l' 


(32b) 


Q  +  “j-k  V* 

m .  ,  1 
0-k  p 


r  ,  C'  ^j+(k+l) 

bj+(k-l)  ml 

j-k  p 


Aside  from  the  trigonometric  terms,  sin  2a  and  cos  2a,  the  coefficients  of 

2 

Eqs.  (32a)  and  (32b)  are  constants.  In  view  of  the  magnitude  of  od^  ,  which 
multiplies  these  functions,  a  series  solution  in  terms  of  a  small  parameter 
will  be  assumed.  Thus 

dAa 


Aa(t,e)  =  Aa  _  + 


e=o  V  de 


€  +  ... 


e=o 


^j*k^e-o  +  (  ~“3e  )e_o  €  +  *** 


(33a) 

(33b) 


where 


CD 


2  q  ? 

ft  +  2  ft  CD  +  7T  (JD 

c  2  c 


2 

CD 

_  C 

ft2 


With  e  equal  to  zero  Eq.  (32b)  is  homogeneous,  and  may  be  solved  in  terms 
of  trigonometric  functions.  After  the  £'s  with  e  zero  have  been  determined, 
Aae_o  can  then  be  found  from  Eq.  (32a).  The  perturbation  differential 
equations  are  obtained  by  differentiating  Eqs.  (32a)  and  (32b)  with  respect 
to  e,  and  then  setting  e  equal  to  zero. 
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3  2  .  ^j-1  "  ?j+l7=o 

AaU  +  2  »  Sln  2ae=o  =  £1  - 

o  p 


I  0  s  cos  2a 
o 


(34a) 


ft'  -  t'  )  2 

^j-1  ^j+3/e=o  Q 


2  1 


^j*k^e=o  +  ^j+k^e=o 


2  2  (  20  +  VkVk 

"  fl  +  Q  (  '■  XD.  ,  1 

J-k  p 


(54b) 


(f  -  ) 
^j+k;€=o 


3  2  3  2 

^  a  cos  2a  +  ^  si  cos  2a 
2i  €  =0  2  €  =o 

o  o 


2Q  4  -.i-kvk 

m .  ,  ,1 
J-fc  p 


+  I  n  2  cos  2a 
o 


j*k 


2  1  ^j-1  "  ^j+l^e=o 
P 


J+l* 


_.Q.  ,  t  +  Q  +  ^j-^j-k  . 

mj-k  Xp  0+(k+l)e=o  m^_k  1  sj*(k-l)e=o 


+  SI 


(  £  I  _  £  I  )  +  Q  f  T 

2  1  '^j-1  ^ j+l  t=o  m.  .  1  ^j*(k+l)e=o 

p  j-k  p  °  ' 


e  +  rn  TJ 

i  _ J j  y  t 

“j-k1?  ;^(k-l)c=o 


where  the  prime  denotes  the  operation  ^  . 

An  examination  of  Eqs.  (34a)  and  (34b)  reveals  a  set  of  non-homogeneous 
equations,  with  the  forcing  functions  arising  from  the  e  equal  zero  solutions. 
Higher-order  perturbation  equations  may  be  obtained  by  continued  differentia¬ 
tion  of  Eqs.  (34a)  and  (34b)  before  setting  e  equal  to  zero. 
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Case  II 

It  is  assumed  that  the  wire  is  nearly  aligned  with  the  local  vertical. 
Under  these  conditions  both  a  and  a  as  well  as  the  (;rs  are  small.  With 
these  approximations  Eqs.  (29a)  and  (29b)  may  be  written  as  follows: 


2  ^i-1  ‘  ^j+1^ 

a  +  3o,c  a  =  . -y'  . 


a 


(35a) 


+  ^  j^k 


.  i m * -i- v» 


m 


i-k^p 


where 


a 


iTk  q  +  m 

(r  -  r  )  + - si 


Q  +  ID<-kVk 

1  u*(k-l) 


J-k  p 


+  Q  ^ J*(k+l) 
mj-k  XP 


(35b) 


a 


n 


+  2  Q  cu  +  3gd 
c  c 


+  2  o> 


For  this  case,  the  set  of  equations  represented  by  Eq.  (35b)  may  be 
solved  immediately  in  terms  of  sinusoidal  functions.  Once  the  £'s  are 
known  Eq.  (35a)  may  be  solved  for  a. 

To  illustrate  the  procedure  described  in  the  preceding  paragraphs, 
a  specific  case  will  now  be  examined.  Let  the  number  of  particles,  n,  be 
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5,  then  it  follows  that  j  is  3  and  the  index  k  is  either  zero  or  one. 
Upon  consideration  of  Eqs.  (29a)  and  (29h)  and  Eqs.  (27a)  and  (27b),  the 
following  expressions  may  be  deduced: 


a  +  3coc  sin  a  cos  a  = 


(C2  -  £4) 

21  si 


(36a) 


^2  +  ^2 


4m  +  m 

-  S2  +  I  S1  +  (— >  S1 


,  . .  fgmx :  m2  - % .  it 

6-1  '  ^3  '  mg  2 

(36b) 


S  +  b 


r  4m,  +  2m 

-  s  +  (— i - 2)  8. 

p  v  hu  1 


2m  +  m 

■  S1  <  V  >^2  +  ^ 


(36c) 


^  +  ^ 


4m  +  m 

-  Ss  +  S  si  +  (— 5T~)  si 


=  S. 


2m  +  m 

c5  <— 


IT?,  .k 


nr 


2 

(36d) 


An  inspection  of  Eqs.  (36a-d)  suggests  that  the  following  transformations 
would  be  useful: 


t2  -  Eu  -  *  <57») 

Eg  +  C4  -  y  -  (57b) 


Thus 


,  2  .  1 
a  +  3oo  sin  a  cos  a  =  n  , 
c  cl 


X 


(38a) 
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x  +  x 


4m  +  m 

si  -  s2  +  ( — 5T — >  si 


=  0 


(38b) 


y  +  y 


4m  +  mp 

2S1  -  S2  +  (  m  "  >  S1 


2m  +  m  -  m 

■  £8i (  m£  . ?)  b 


(58c) 


V53 


4m.  +  2m_ 

-  *2  +  S1 


2m  +  m 

-  sx  (  7  s)  y 


m. 


(38d) 


We  are  now  in  a  position  to  obtain  the  solutions  of  Eqs.  (38a-d)  for 
the  two  cases  under  consideration. 

Case  I 

The  first  step  is  to  express  x,  y  and  £  in  a  Taylor  series  in  e. 


Act  =  Aae_0  + 


y  de  /e=o 


e  +  . . . 


dx 

e=o  \  SI  i  6=o 


(39a) 

(39b) 


y  -  ye=o  +  1ST  Js=o 


€  +  .  .  . 


(39c) 


h  =  (^3)6=0  +  \^r)^C 


e  +  . . . 


(39d) 


With  e  equal  to  zero  Eqs.  (38b-d)  may  be  solved  immediately. 


x  =  A,  cos  (Jo  t  +  A  sin  oo  t 
e=o  1  x  2  x 


(40a) 
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The  A's,  B's  and  C's  are  functions  of  the  initial  values  of  the  ^'s  and 

•  t 

the  £'s.  For  convenience  it  will  be  assumed  that  all  of  the  £'s  are 
initially  zero.  Under  these  conditions, 

A1  =  x(o)  =  C2(o)  -  ik(o) 
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+  £5(o)  q 


2 

o 


-  S(o)  °o 


+  2m ^  -  2m ^ 


1  2 m  +  m 

!)  +  y(°)  n0  ( — e—  > 

-J  n 


p  2m  +  m 

-  y(°)  n0  ( - m - ) 


of  Eqs.  (38)  are  oscillatory  in 
t  expressions  for  co^  and 
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indicates  that  for  any  mass  distribution  that  is  physically  possible,  the 
e  equal  zero  solutions  are  always  composed  of  simple  sinusoids. 

With  x  known  it  is  now  possible  to  find  A2  __  from  Eq.  (38a). 

6=0  t-O 


to. 


da 


e=o 


21 


A.,  sin  cd  t 
1 _ x_ 

0) 


A„  COS  CD  t 
2 _ x_ 

CD 


A  n2 

A2  no 

21  CD 
P  X 


(hi) 


We  can  now  check  the  validity  of  the  assumption  that  da  is  small 
and  of  the  order  of  x.  With  A^  zero 


Aa 


€=0 


ft  x(o)  sin  cd  t 
o  _ x 

^ml  +  m2 


21 


m„ 


2 

Due  to  x(o),  da  and  x  da  are  second-order  terms.  For  a  wire 
6=0  e=o  €=0 

2000  ft  long  and  with  x(o)  equal  to  1  cm,  the  magnitude  of  Aa£=o  is 

approximately  3  x  10  ^  radians  per  second.  Thus  the  retention  of  the 

cd2  terms,  in  s  and  s0,  as  compared  to  the  Q da  is  justified  only  if  x(o) 
c  12 

-3 

is  on  the  order  of  10  cm. 

The  first-order  perturbation  equations  are  found  by  differentiating 
Eqs.  (38a-d)  with  respect  to  e,  and  then  letting  e  go  to  zero. 


da'  +  \  ft2  sin  2a 

6=0  2  O  6=0 


6=0 
2  1 


3  2 

(|  fl0  cos  2 Qe=0)  +  2Y 


6=0 


(42a) 


X'  +  X' 
6=0  6=0 


CD 


=  -  X 


6=0 


3ft  cos  2a  + 
o  e=o 


4®,  +  m2  x  P 

( — - -)  5  cos  2a 

m  '2  0  €=o 


(42b) 
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^  2  2 

y '  +  y’  (ft  +0))  = 

Je=o  Je=o  '  o  x' 


ye=o 


"q  p  H  +  m2  ^  2 

I  S>o  c°s  2a£=o  +  (  „— )  I  ft0  cos  2ae=o 


2m  +  m  -  m 

+  31?  cos  a a£=o  ( - — - 2)  (Cj)£=0 


(42c) 


2m  +  m  -  m 


r  „  km,  +  2m„ 

~ 

+  ^p€=0 

^2/1  „2 

-  ft  +  ( - )  ri 

o  m,  'o 

5 

^3^€=o 

l  *o  C°S  ^€=0 


V  +  *  p 

+  < - 5—)  I  C°S  ^6=0 


,  o  2m  +  m 

+  i  ft0  cos  ^=0  (  m:~-)y€=o 


+  m„ 


m. 


(42d) 


Equations  (42a-d)  are  linear,  but  non-homogeneous .  Since  all  of  the 
initial  conditions  have  been  satisfied  by  the  e  zero  solutions,  the  initial 
values  of  displacement  and  rate -of -change  of  displacement  must  be  zero  for 
all  of  the  perturbation  variables. 

Upon  substitution  of  Eqs.  (40a-c)  into  Eqs.  (42a-d)  the  solutions  of 
Eqs.  (U2b-d)  may  be  found 


x 


e=o 


x(o)  (5  +  |  a?) 

8ft  (od  +  ft)(o3  -  ft) 

X  X 


2ft  COS  03  t 
X 


(o3  -  ft)  COS  (2ft  +  03  )t 

X  X 
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-  (a>x  +  a)  cos  (2a  -  <^x)t 

y'  =  B’  cos  a>nt  +  BI  cos  as„t  +  K,  cos  (2a  -  as  )t 

J  e=o  1  1  5  2  1  x 

+  Kg  cos  (20,  +  u>^)t  +  K^  cos  (2a  -  a)g)t  +  K^  cos  (2a  +  0Sg)t 

(f  l)_  =  C’  cos  0Dnt  +  Ci  cos  0)ot  +  Kc  cos  (2a  -  as,  )t 

Vb3/€=o  1  13  25  1 

+  Kg  cos  (2a  +  as1)t  +  cos  (2a  -  cUg)t  +  Kg  cos  (2a  +  a>g)t 


(43a) 

(43b) 


(43c) 


where 


2  2 
as2  -  asx 


3n2  (Cx  +  C5)  -  +  B2^  "  “2  ^K1  +  K2  +  K3  +  K4^ 


+  Kx(2a  -  a^)2  +  K2(2a  +  as^2  +  K^(2a  -  cd2)2  +  K^(2a  +  a>2)2 


■R 1  — 

3  2  2 
mi  ’  “2 


3n2  (C1  +  C  )  -  n1(B1  +  B  )  -  as2  (K±  +  Kg  +  K  +  K^) 


+  Kx(2a  -  as.^2  +  Kg  (2a  +  a^)2  +  K^(2a  -  aSg)2  +  K^(2a  +  as2)2 


r '  =  — - - 

1  2  2, 

“2  ’  “i  L 


|  (B1  +  B?)  -  n3(C1  +  C?)  -  as2  (K$  +  Kg  +  K?  +  Kg) 


+  K  (2ft  ~  as^)  +  Kg  (2a  +  as^)  +  K^  (2a  ~  0Sg)  +  Kg (2a  +  aSg) 
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C3  2  2 
“l  “  “2 


I  n4  ^Bi  +  "  n3^Ci  +  C "  “l  ^K5  +  k6  +  K7  +  K8^ 


+  K  (2a  -  o^)2  +  Kg(2ft  +  o^)2  +  K?(2ft  -  cd2)2  +  Kg(2ft  +  a>2)‘ 


3n2  n4  ~  2n2  n3  ^1  +  ^2.  "  nlBl^ 

4m  +  2m  -  m  P 

(  1  „  2  5to2  -  (a  +  »  )2 

m  °  1 

L  p  -l 

8ft  (XV  -  £D1)  (2ft  “  0^  +  GQ2)(20  -  -  o>2) 


5n2  ni^  ”  ^n2  +  ^5n2  0-^  - 


K. 


4m  +  2m  -  m  2 

(  1  „  2 - ^  -  (a  +  cd1)2 


m. 


8ft  (ft  +  a^)  (2ft  +  oi1  +  oj2 )  (2ft  +  od1  -  (X>2 ) 


3n2  n4  B3  ‘  2n2  n3  C3  +  (3n2  C3  '  niB3^ 


4m  +  2m  -  m  p  2 

(  1  ■  2  -  (a  -  “ft 


8ft  (ft  -  CD2)(2ft  -  CJJ2  +  0>1)(2ft  -  a>2  -  cox) 


K. 


4m,  +  2m„  -  m„  _  _ 

5n2  n4  "  ^n2  n3  ^3  +  ^5n2  ”  nl  B3^ 

(  1  m  g — -  (a  +  «fe> 

= _ 2 _ I - J 

(ft  +  co2)  (2ft  +  cd2  +  co1)(  2ft  +  co2  -  a>x) 


3n2  P4  C1  ~  ni  n4  B1  +  ^1  \  B1  *  n3  Cl^ 


2  2  2 
ftQ  +  <*>x  “  (2ft  -  cdx) 


8ft  (ft  -  CD1)(2ft  -  +  a>2)(2ft  -  -  cd2) 


3n 


K, 


2  n4  C1  nl  n4  B1  +  ^2  n4  E1  "  n3  Cl^ 


2  2  2 
ftQ  +  a>x  -  (2  ft  +  oo^) 


8ft  (ft  +  ax^)(2ft  +  ax^  +  co2)  (2ft  +  -civ,) 
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?n2  C3  -  nx  B;  +  (g  ^  B?  -  C?) 


2  2  2 
ftQ  +  cox  -  (2ft  -  cOg) 


8ft(ft  -  C0g)(2ft  -  o)g  +  cu1)(2fi  -  ciig  -  oi1) 


Kr 


3ng  *4  c3  -  D1  n4  B3  +  (J  n4  B3  -  n3  C^) 


2  2  2 
00  +  0>x  -  (2ft  +  03g) 


8ft  (ft  +0Dg)(2ft  +  a> g  +a>1)(2ft  +  a>g  -  ) 


2  mi  + 

“i  -  *  (-hr1)  * 


2m  +  m  -  m 

n  «  ft2  (— ± - 2 - 2) 

2  o  m.  ' 


,  2  ^“l  +  2m2  + 

n3  =  2  Ao  ^  ml 


)  ; 


.2  ,2ml  +  m2, 

n4  "  fto  ( - i— > 


Finally,  the  perturbation  spin  rate  may  be  found  by  introducing  Eq.  (43a) 
into  Eq.  (42a). 


A“£=o 


3  o 
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1)  ’  r  - p~  sin  20 1 

1  to  SI 
P  x 
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An  examination  of  the  first-order  perturbation  solutions  indicates 
that  a  resonance  condition  exists  if  any  one  of  the  following  conditions 
is  met: 


05^  -0=0  (a) 

-  O  =  0  (b) 

CD  -  0  =  0  (c) 

20  -  03  +  o>2  =  0  (d) 

2Q  -  cd^  -  co2  =  0  (e) 

20  -  o>2  +  =  0  (f ) 


(45) 


Case  II 

When  the  five-particle  wire  is  aligned  with  the  vertical,  Eqs.  (35a) 
and  (35b)  may  be  written  as  follows: 


^  3co 
...  2  >UJc 

a  +  3o3c  a  =  x 


(46a) 


x  +  >c  ( 


2  +  ^ 


)  x  =  0 


(46b) 


4m  +  3m 

y  (  1  2 


- )  y  =  6cn  (— - - - 2)g 

nu  '  *  c  m_  'b3 


(46c) 
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4m 

(— 


+  2m, 


EL 


+  m. 


m. 


^y 


(46d) 


Equations  (46a-d)  are  linear  with  constant  coefficients,  and  thus  may  be 
solved  in  terms  of  exponential  functions . 


/ 

X  =  x(o)  COS  l  3( - 


4m,  +  2m, 


m. 


Ot 


(^7a) 


a  = 


a{o)  +  — hgi 


21 


cos 


vT  coct 


X 


’o) 


~ 4m, 

2 1  /  1\ 

p  (_£T) 
2 


cos  CO 


4m,  +  2m 
3( — =:  2 


m. 


)t 


(47b) 


y  =  B^  cos  ayt  +  cos  co^t 


(47c) 


C.  =  C,  cos  cd  t  +  C,  cos  a3„t 

b5  i  y  3  £ 


(47a) 


where  the  B's  and  C’s  are  functions  of  y(o)  and  ^(o)  and  where 


2  2 

co  =  3u> 

y  c 


4m^  +  2m2  41^  +  3m2 


2m. 


2m, 


4m  +  2m  2  4m  +  3m  2  4m  +  2m 

<  ■  V  -)  +  < - a^)  -  <-  ay3-  )(m2  +  ^ 
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An  examination  of  Eqs.  (47  s-d)  indicates  that  for  small  disturbances, 

the  wire  oscillates  about  the  vertical.  Without  damping,  these  oscillations 

continue  indefinitely,  but  the  vertical  configuration  of  the  wire  is 
* 

completely  stable. 

*If  the  orbit  has  an  eccentricity  much  greater  than  zero,  the  above 
conclusions  are  not  valid. 
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IV.  DISCUSSION  OF  THE  RESULTS 

The  results  of  the  previous  section  indicate  that,  for  the  spinning 
wire,  parametric  resonance  can  occur.  In  the  example  analyzed,  only 
three  normal  modes  of  vibration  are  possible.  (See  Fig.  5)  The  center 
particle  can  oscillate  at  only  two  different  frequencies,  oo^  and  cd^.  The 
remaining  particles  can  oscillate  at  three  different  frequencies,  oj  , 

■X- 

and  Thus,  because  of  the  finite  number  of  frequencies,  resonance 

will  occur  only  if  one  of  the  six  possible  combinations  of  the  of  s  with  0 
happens  to  be  exactly  zero.  However  as  the  number  of  particles  increases, 
the  number  of  conditions  under  which  resonance  takes  place  also  increases. 
Thus  as  n  becomes  infinite,  resonance  would  surely  occur,  for  the  mathe¬ 
matical  wire  under  study.  The  initial  rotational  kinetic  energy  would 
eventually  be  transferred  to  the  various  vibrational  modes  of  the  wire. 


V 


Fig. 5 — The  normal  vibrational  modes  for  five  particles 


All  of  the  particles  rotate  with  the  tj  -  £  axes  at  the  angular  rate, 
(a  +  a>c). 
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It  is  doubtful  that  a  real  wire  would  be  rotationally  unstable  in 
the  sense  just  described.  In  order  to  see  this  more  clearly,  let  us 
review  the  physical  and  mathematical  approximations  that  have  been  adopted 
in  the  development  of  the  equations  of  motion. 

Physical  Approximations: 

1.  Assumption  of  a  spherical  gravitational  potential  function. 

2.  The  neglect  of  internal  restoring  and  dissipation  forces,  i.e., 
spring  constant  and  damping  terms. 

3.  The  neglect  of  all  external  forces  with  the  exception  of  number  1. 
Other  forces  which  act  on  the  wire  include  those  due  to  inter¬ 
action  with  the  earth’s  magnetic  field  and  electrostatic  forces. 

4.  Restriction  of  the  motion  to  a  plane. 

Mathematical  Approximations: 

2 

1.  The  neglect  of  terms  such  as  ^  etc. 

2.  The  neglect  of  terms  of  second  degree  and  higher  in  the  expansion 
of  r^« 

Of  the  physical  approximations,  number  2  is  probably  the  most  signifi¬ 
cant  from  the  viewpoint  of  stability.  In  particular,  the  addition  of  dis¬ 
sipative  internal  forces  would  modify  the  results  of  the  preceding  section* 
For  a  continuous  wire,  such  damping  forces  arise  due  to  internal  friction 
between  intergranular  boundaries.  An  exact  analysis  of  the  problem  is 
difficult  due  to  the  fact  that  the  damping  terms  can  be  functions  of  both 
the  amplitude  and  the  rate  of  change  of  the  £'s.  As  a  consequence  it  is 

*  ( 5 ) 

even  possible  that  the  damping  may  be  a  destabilizing  effect. 

*In  essence,  the  system  might  be  self-excited  due  to  an  effective 
negative  damping  coefficient.  The  energy  is  supplied  by  the  over-all 
rotational  motion. 
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In  the  Appendix,  the  effect  of  viscous  damping  on  the  motion  of  the 

x  normal  mode  is  considered.  It  is  shown  that  very  small  damping  ratios, 

on  the  order  of  0.05,  are  more  than  sufficient  to  eliminate  any  divergence 

problem.  Furthermore,  as  long  as  the  initial  disturbances  are  such  that 
2 

the  neglect  of  £  terms  is  valid,  steady  state  conditions  with  essentially 
no  transverse  deflections  are  reached  within  a  few  minutes. 

The  neglect  of  spring-restoring  forces  is  a  reasonable  approximation, 
since,  for  a  flexible  wire,  they  are  very  small  compared  with  the  tension 
restoring  forces.  If  they  were  included  in  the  analysis,  the  values  of  the 
undamped  natural  frequencies,  o^,  a etc.,  would  be  slightly  increased. 

Another  physical  approximation  which  might  be  important  is  number  3. 

A  conducting  wire  rotating  in  the  earth’s  magnetic  field  would  be  subjected 
to  sinusoidally  varying  forces.  An  external  forcing  function  of  this  type 
would  be  of  importance  if  a  resonance  condition  occurs. 

Of  the  mathematical  approximations,  number  1  is  probably  the  most 
important.  The  coupling  between  the  vibrational  and  rotational  degrees  of 
freedom  is  nonlinear  in  nature.  Thus  the  linearization  procedure  prior  to 
the  perturbation  analysis  eliminates  the  forced  solution  due  to  terms  such 
as  x  a a.  It  does  not  appear  that  the  gross  characteristics  of  the 
vibrational  motion  of  the  wire  are  altered  by  the  neglect  of  such  terms. 
However,  when  the  forced  motion  of  the  wire  with  damping  is  considered, 
the  energy  dissipation  of  the  system  is  strongly  influenced  by  the  non¬ 
linear  coupling  terms.  As  a  consequence,  estimates  of  the  lifetime  of  the 
spinning  wire  cannot  be  obtained  from  the  linear  model  that  has  been 


examined. 
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The  primary  causes  of  forced  motion  are  the  interaction  with  the 
earth's  magnetic  fields  mentioned  above,  and  non-circular  orbits. 

Other  external  forcing  functions  would  appear  on  the  right-hand  side  of 
the  equations  of  motion  if  higher-order  terms  are  retained  in  the  expansion 
of  r^,  and  if  the  earth's  oblateness  were  considered. 

The  introduction  of  damping  eliminates  the  possibility  of  either 
parametric  or  externally  forced  resonance.  However,  the  continuous 
excitation  of  the  wire  due  to  such  second-order  effects  as  an  oblate  earth, 
a  non-circular  orbit,  rotation  of  the  wire  out  of  the  orbital  plane,  etc., 
will  cause  a  gradual  dissipation  of  the  combined  rotational  and  vibrational 
kinetic  energy.  Thus  the  spin  rate  of  the  wire  will  gradually  decrease. 
When  the  tension  forces  due  to  spin  are  of  the  same  order  of  magnitude  as 
the  gradient  tension  forces,  the  straight  line  configuration  of  the  wire 
will  no  longer  be  stable.  • 
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Appendix 

THE  INFLUENCE  OF  VISCOUS  DAMPING 

As  indicated  in  Section  IV  the  internal  dissipative  forces  of 
a  wire  cannot  be  represented  by  a  simple  viscous  model.  However,  it 
might  be  possible  to  coat  the  wire  with  a  substance  which  would  impart  the 
desired  damping  characteristics.  In  order  to  examine  the  influence  of 
dissipative  forces  on  the  motion  of  the  wire,  it  is  assumed  that  a  viscous 
damping  process  is  applicable.  From  Fig.  2 

Fai~  si-i,i  '  pi,i+i  “  +  4+iJ  (48) 

If  it  is  assumed  that  the  damping  is  uniform  per  particle,  then  for  the 
five-particle  case 


F, 

dl 

-  2  - 

(49a) 

% 

-  a(q  +  -  2(g) 

(49b) 

s 

-  d((2  +  (4  -  2(j) 

(49c) 

F, 

d4 

-  d(c5  +  c5  -  2C4) 

(49d) 

Fd 

e: 

-  -  (5) 

(49e) 

'5 


Utilizing  the  center  of  mass  relationships  to  eliminate  and  t,  ,  the 
equations  of  motion  are 
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“  +  a  “c  sln  aa  =  £1  +  m2y  +  21;  x  S1 
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x  +  d 
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x  +  x 


,*“1  +  a 
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h  +  2  m  S  +  S 
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4m1  +  2m 
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2  v  m_  '  j 
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The  solutions  of  Eqs.  (50a-d)  can  be  found  by  the  methods  employed 
in  Section  III.  To  illustrate  the  effect  of  damping,  however,  we  will 
only  consider  Eqs.  (50a)  and  (50b).  With  x(o)  =  0,  the  €=o  solution  of 
Eq.  (50b)  is 


-£co  t  to 

x  =  x(o)e  x  (cos  cot  +  |  —  sin  cot) 
e=o  ^ 


(51) 


where 


co 
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d  ^  (2mx  +  m2)  - 

m1m2(2m1  +  mJJ 


If  it  is  assumed  that  the  damping  ratio,  is  much  less  than  one. 


^  ft2  x(o) 

aA  ~  _2 -  g  +  (---4?  sin  tot  -  i  cos  cot)  e 

e=o  1  to  2  to 

-Lp  x  L  x 


-£co  t 
x 


dx(o) 


V7.^ - \  (e  *  cos  tot  -  1) 

l(2m1  +  m2)  v 


The  x  perturbation  equation  of  motion  may  be  found  by  differentiating 
Eq.  (50b)  with  respect  to  e  and  then  setting  e  equal  to  zero. 


••  ’2  2  5  2  "*»* 

x’=Q  +  2|o>x  x*=o  +  tox  x»=o  =  -(5lo  +  |  tox)  x(o)  e  cos  tot  cos  2ftt 


+  -  sin  tot  cos  2">t 

to 


The  solution  of  Eq.  (53)  is 


-£co  t 


x'  =  e 
e=o 


Aj^  cos  cot  +  A2  sin  tot 


x(o)(3a^  +  |  co2)  cos  (2fl  -  co)t 
77  (co  -  fl) 


cos(2fl  +  co)t  |co  sin  (to  +  2Jl)t  |co  sin  (to  -  2ft  )t 

x  _ _ _  +  _ 2S  _ _ 

(to  +  ft)  co  (co  +  ft)  to  (co  -  TT~ 


(54) 
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With  x*  (o)  =  x^=o(o)  =  0,  Eq.  (54)  becomes 


x'  =  e 
e=o 


-Soyt  x(o)(3fl2  +  |  ni2)  r 
Wl  co-fl)  (co+fl) 


2Q,  cos  cot  -  (a>+.ft)  cos  (2fi-co)t 


So) 

+  (co-fi)  cos  (2fi-ko)t  -  2ft  sin  cot 


(55) 


£C0  £C0 

+  (co-ft)  — -  sin  (co  +  2ft)t  -  (co+ft)  — -  sin  (co-2ft)t 
v  '  05  0) 


As  Q  approaches  co,  resonance  occurs.  Thus 


x(o)(3fi2  +  I  <o2)  -£co  t 
v  2  x'  x 


e=o 


16 


co 


cos  5  cot  -  cos  oyt  -  4ot  sin  cot 


£o>x 

+  -  (sin  cot  +  sin  5  cot)  -  4£co  t  cos  cot 


CD 


(56) 


The  total  solution  for  the  normal  vibrational  mode,  x,  with  viscous 
damping  included,  is: 


i  2.  252s 

-|cot  co  co  (3ft  +  ?  co  ) 

/\  x  J  „  x  .  ,  c  '  o  2  x' 
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16  co  ft 
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(57) 
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If  circular  orbits  with  altitudes  of  the  order  of  200  to  2000  miles  are 
employed,  then,  without  damping,  the  secular  term  of  Eq.  (57)  predominates 
after  a  period  of  ten  days  to  two  weeks.  However  even  with  a  value  of  | 
of  the  order  of  0.05,  the  amplitude  of  x  is  l/e  of  x(o)  at  the  end  of  20 
seconds.  The  secular  terms  are  completely  dominated  by  the  exponential 
after,  approximately,  100  seconds. 

Differentiation  with  respect  to  e  of  Eq.  (50a)  yields  the  a 
perturbation  equation. 


'“ho  ■  tt  (3/2  no cos  2rt)  +ri 

P  P 


r  t 

6=0 


(58) 


d  x  ’ 

— - - - - -  -  \  f>2  sin  (2Qt  +  2  to.  ) 

1  (Am  +  2m  )  2  o  €=o/ 

T)  f  k 


.  t  .1 

With  the  substitution  of  x  ,  x*  ,  x  and  Ad  into  Eq.  (58),  to. 

e=o  e=o  £ -O  £=0  e=o 

may  be  found  by  a  direct  integration.  The  resulting  expression  is  quite 
complex,  but  it  can  be  easily  shown  that  all  of  the  terms  quickly  damp  to 
zero  with  one  exception.  This  particular  term  would  exist  even  if  the  body 
was  rigid  since  it  arises  from  the  gravitational  gradient  torque  term  of 
Eq.  (58).  Physically  the  effect  may  be  explained  as  follows:  As  the  wire 
rotates  about  its  center  of  mass,  the  potential  energy  of  the  wire  changes 
periodically.  The  change  in  potential  energy  is  reflected  in  a  change  of 

the  rotational  kinetic  energy  of  the  wire.  Thus  the  spin  rate  also  varies 

* 

periodically.  This  variation  in  the  spin  rate  is  extremely  small,  and  it 

2 

has  no  effect  upon  the  vibrational  motion  until  terms  of  the  order  of  e 

are  retained. 

* 

An  exact  analysis  of  this  effect  requires  that  the  coupling  between 
the  rotational  motion  and  the  orbital  motion  be  included  as  indicated  in 
Ref.  3. 
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